Series of Abelian and Non-Abelian States in C > 1 Fractional Chern Insulators 
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We report the observation of a new series of abelian and non-abelian topological states in fractional 
Chern insulators (FCI). The states appear at bosonic filling v = k/(C + 1) (k, C integers) in a wide 
variety of lattice models, in fractionally filled bands of Chern numbers G > 1 subject to on-site 
Hubbard interactions. We show strong evidence that the k = 1 series is abelian while the k > 1 
series is non-abelian. The energy spectrum at both ground-state filling and upon the addition of 
quasiholes shows a low-lying manifold of states whose total degeneracy and counting matches, at 
the appropriate size, that of the Fractional Quantum Hall (FQH) SU(C) (color) singlet fc-clustered 
states (including Halperin, non-abelian spin singlet(NASS) states and their generalizations). The 
ground-state momenta are correctly predicted by the FQH to FCI lattice folding. However, the 
counting of FCI states also matches that of a spinless FQH series, preventing a clear identification 
just from the energy spectrum. The entanglement spectrum lends support to the identification 
of our states as SU(C) color-singlets but offers new anomalies in the counting for C > 1, possibly 
related to dislocations that call for the development of new counting rules of these topological states. 



PACS numbers: 73.43.-f, 71.10.Fd, 03.65. Vf, 03.65.Ud 

In his seminal paper [T], Haldane introduced the con- 
cept of a Chern insulator (CI). Classified by a non-zero 
Chern number of the occupied bands, this band insula- 
tor exhibits an integer Hall conductance akin to the inte- 
ger quantum Hall effect but in the absence of an overall 
magnetic field. Recently, the role of strong interactions 
in CI has received strong attention [2T-H3] focused on re- 
alizing the FQH in flat bands akin to Landau levels but 
in zero field. Most of these studies have focused on CI 
with fractionally filled bands of Chern number C unity. 
These systems have been dubbed Fractional Chern Insu- 
lators (FCI). Strong numerical and analytical evidences 
of Laughlin-like 13] states[2 6,15 or more generally com- 
posite fermion[16 states have been obtained in CIs in the 
presence of 2-body interaction [17] . With (k + l)-body 
interactions, phases similar to the Moore-Read 18J (MR) 
or Read-RezavifTD] (RR) states emerge [31 |5"II2"U]. 

A nice feature of CIs absent in the QH effect is the 
possibility to obtain a band with a higher Chern number, 
without introducing an additional degree of freedom (like 
the spin in the FQH). Several simple one-body models of 
C > 1 CIs with an almost flat band now exist [2"TM25| . 
Ref. [35] proposed [2DJ , based on the Wannier approach 
[T5"ll2"Tll2"8] . that non-interacting CIs with C > 1 are iden- 
tical to the the quantum Hall effect with particles carry- 
ing an extra SU(C) internal degree of freedom. It has 
then been assumed that this property should also hold 
true in the strong interaction regime. Moreover, several 
recent numerical works [531 UHl EDJ have reported the ob- 
servation of FQH-like phases at filling factor v = for 
bosonic systems and v = 2 c+i f° r fermionic systems, and 
presented evidence (based only on the energy spectrum) 
that they are abelian. While the possibility that these 
states have an internal hidden structure (such as some 
SU(C) or Z c degree of freedom) was raised in 23, 29, 30J, 



no evidence for this scenario was provided, and the ob- 
served energy spectra could equally be well-explained by 
spinless states. 

In this paper, we pose the natural question of whether 
generalizations of the MR or RR states can emerge in FCI 
when C > 1. In this letter, we present strong indications 
that such phases are realized in FCI at filling factor v — 
P^-j- when a fractionally filled CI band is subject to a 
(k + l)-body interaction. If the analogy between the FCI 
and the SU{C)FQH is valid, we expect that this series 
to be related to the Halperin state[3T] (for k = 1), the 
non-abelian spin singlet [351 (NASS) (for C = 2) states 
and their generalizations [33] (for higher k values). 

The signatures of the FQH-like phases in FCIs ap- 
pear in different manners [JJ, the most simple of which 
are the groundstate (quasi-)degeneracy, the counting of 
the number of quasihole states and their momentum 
quantum numbers whose values can be obtained from 
the FQH-FCI mapping described in [3J. Unfortunately, 
a Charge Density Wave (CDW) would exhibit similar 
counting[3I] in its energy spectra. The particle entan- 
glement spectrum[35l [36] is a better way to identify sig- 
natures of FQH phases [37], and also to distinguish be- 
tween CDW and FCI. The previous studies[H HH [3DJ 
of C > 1 FCIs have mostly focused on the matching be- 
tween counting of groundstate or quasihole manifold de- 
generacy of the FCI and FQH. Even excluding the pos- 
sibility of CDW phases, as we clearly show this is not 
enough to prove that the physics of C > 1 FCIs is re- 
lated to the spinful FQH effect: we show that the count- 
ing of any of the previously observed states [53J (3DJ [3DJ 
or our new non-abelian series can be deduced from a 
spinless generalized Pauli exclusion principle [38j [39]. We 
show that the entanglement spectrum is able to rule 
out the option of a spinless FQH states while providing 
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indication (some yet un-understood properties disallow 
a clear proof) of colorful FQH states. More recently, 
approaches [El [40] [41] , some based on the gauge- fixed 
unitary |15j Wannier basis [2"Tj were used to compute over- 
laps and to adiabatically continue between FCI and FQH 
states in C = 1 FCIs, but this approach is not yes appli- 
cable in the current C > 1 case |26| . 

In this letter, we mostly focus on the slightly simplified 
pyrochlore lattice model introduced in (33], consisting of 
N Kagome layers coupled through intermediate trian- 
gular layers. For completeness, in the notation of [2~4] . 
we use parameters ti = 1, t± = —1.03, Ai = 0.83 and 
^2 = A2 = (i.e. we have discarded the second nearest 
neighbor hopping). This model gives a flat band with 
Chern number C = N . We also investigate other CI 
models: the triangular lattice model [53] with C = 2, the 
two orbitals on a triangular lattice [23] with C = 3 and 
the N orbitals on a square lattice with C = N [25] . 

We consider only interacting bosons. For the py- 
rochlore and the triangular lattice models, we use a 
(k + l)-body Hubbard interaction H- lnt ,k — UJ2i '■ Pi +1 '■ 
where :: denotes the normal ordering and the sum runs 
over all the sites. For N-Orbital on square lattice and two 
orbitals on triangular lattice models, we use an isotropic 
on site interaction H[ nty k = Uj^i '■ (J2jPij) k+1 '■ where 
the first sum runs over all the sites whereas the second 
sum runs over the different orbitals on the same site. 



Following Ref. [37] , we remove the spurious effect of band 
dispersion and band mixing by using only the Hilbert 
space of the i-th band and neglecting its kinetic energy 
(standard flat band procedure) . We assume that we have 
infinite band gaps and the i — 1 first bands are filled and 
inert. 

As very recently reported in |23j for Chern number 
C = 2 and in j^H] for C = N, there is clear evidence of 
Abelian phases at filling v = 1/(C+ 1) for bosons with 
#int,fc- The energy spectra for C — 2,3 are shown in 
Fig. [T^,) and Fig. [TJ;) respectively. For both cases and for 
all numerically accessible numbers of particles, we find 
(C + l)-fold quasi-degenerate groundstates. The ground- 
states momenta can be predicted from the FQH to FCI 
mapping |34| (which technically works only for C = 1 
but which we empirically find to work here as well) from 
the partition obeying the (1,C + l)i (spinless) general- 
ized Pauli principle (we denote (k,r)c the (fc,r) gener- 
alized Pauli principle with SU(C) symmetry to be de- 
fined below in Eq. [3]) . In [25] , a mismatch between the 
FQH-FCI mapping and the ground-state momenta was 
reported. This is due to the use of fermionic momenta 
of partitions instead of needed bosonic ones. The total 
quasihole counting also matches the total number of par- 
titions obeying the (1, C+l)i generalized Pauli principle, 
suggesting an abelian state. 
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FIG. 1. (a): Low energy spectra on the pyrochlore lattice with C = 2 and two body interaction for N e = 6, 7, 8, 10 bosons at 
v = 3. (b): Low energy spectra on the pyrochlore lattice with C — 2 and three body interaction for 7V e = 6, 8, 10, 12 bosons at 
v = 2/3, we observe an almost sixfold degenerate ground state only for an even number of particles, (c): Low energy spectra 
on the pyrochlore lattice with C — 3 and two body interaction for N e ~ 6, 7, 8, 9 bosons at v = |. (d): Low energy spectra on 
the pyrochlore lattice with C = 3 and three body interaction for N e = 6, 8, 10 bosons at v — 1/2. (e): Low energy spectra on 
the pyrochlore lattice with C = 3 and three body interaction for N e — 7,9 bosons at v = 1/2. We observe an almost tenfold 
degenerate ground state only for even number of particles. We only show the lowest energy per momentum sectors in addition 
to the degenerate groundstate. For C — 3 three body, only the N e even sector exhibits a full gap between the degenerate 
groundstatate manifold and the excited states, a signature of pairing. 



Using the three-body interaction (k — 2), we find a 
six- fold quasi-degenerate groundstate at filling v = 2/3 
for C — 2 and a ten-fold quasi-degenerate groundstate at 
filling v = 1/2 for C — 3 and even numbers of particles. 
The energy spectrum for C = 2 is shown on Fig. [T]d) 



while the one for C — 3 is shown on Fig. [TJi). The total 
number of groundstates and quasihole states are com- 
patible with the (2, C + l)i generalized Pauli principle 
[55] , which suggests a paired, nonabelian state. Another 
signature of the paired state comes from the even-odd 
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particle number aliasing for C = 3. While for C = 2 
the groundstate can only be realized with N e even due 
to the filling factor, this constraint does not apply for 
C = 3 for which we observe a gap separating the almost 
tenfold degenerate groundstate manifold from the higher 
energy excitation only for N e even a clear signature of 
pairing. 

For a general (k + l)-body interaction, we expect to 
find a topological phase at filling v = p^j, for which 
the groundstates and quasiholes countings are given by 
the (k,C + l)i generalized Pauli principle. However, in 
every model investigated so far, we have only found weak 
or conflicting evidences (such as corrected groundstate 
degeneracy but no entanglement gap) of these phases for 
k > 2. One should keep in mind that FCI states are 
highly dependent on the underlying tight binding model 
and interaction range. Thus the absence of the full series 
is most probably related to the peculiarities of the models 
used than due to a fundamental reason. Among the other 
models we have studieds, we also observe clear signatures 
for the following states: k = 2, C = 2 in the triangular 
lattice model[33], fc = 1,(7 = 3 in the two orbitals on 
a triangular lattice[3S] and k = 1,(7 = 4, 5 in the N 
orbitals on a square lattice [25], 

For the cases studied, while the total number of 
groundstates and quasihole states is compatible with a 
(fc, C+l)i generalized Pauli principle, differences appear 
in the observed counting per momenta sector: it does not 
match in all cases the counting obtained using the FQH 
to FCI mapping^] developed for (7 = 1 systems. In par- 
ticular, when the Pauli principle applied is an inherently 
"fermionic" one when you consider FQH (as (1, 3) for ex- 
ample), the FQH counting you obtain for bosons in the 
FQH reduced BZ does not have the C4 symmetry. This 
makes the FQH to FCI mapping fail to predict the count- 
ing per momentum sector. This shows the need for the 
development of a new folding mapping for the C > 1 
systems [2"8] , 

While the total counting of degenerate groundstates 
and multiplcts matches a spinless FQH counting, this 
is not evidence that the state does not have some in- 
ternal SU(C) symmetry. To see this, we analyze can- 
didate FQH systems with SU(C) internal symmetry, of 
which the generalized [m; n]-Halperin wavefunctions are 
the main abelian canditates. They are given by: 
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where = Y\ti Ylk i<h ( z % ~ 4?)™ is the prod- 

uct of a Laughlin state for each component and ^V^? = 

I\f<j 11^=1 Ukj=i( z S ~ 4?)" accounts for correlations 
between components. Here, z^ is the complex position 
of the fc-th particle of component i. The exponents m and 
n characterize the strength of the intra and inter compo- 



nent correlations respectively. The total filling factor is 
' and the ground state degeneracy on 



l)n). These states 
1. In our case, we 
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the torus is d — (m — n) ^ 1 (m + (C - 
are Si7((7)-singlets when n = m 
consider m = 2,n = 1. 

From the states in Eq. [I] one can build a series of 
non Abelian spin-singlet states[32 at filling i^fqh = 
c(m-i)+i • This is done by dividing the particles into 
k groups, writing a Halperin [m; m — 1] state for each 
group and then symmetrizing over the different groups. 
This procedure leads to: 
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where S is the symmetrization operator. For m = 2, n = 
1, these are the NASS states introduced by Ref. [32] for 
C = 2 and generalized for C = 4 in Refs. g2J S3J . 

The counting of excitations of the series above has been 
determined through generalized Pauli principles (321 144j . 
All spinless fermionic (bosonic) many-body wave func- 
tions of N e particles can be expressed as linear combina- 
tions of Fock states in the occupancy basis of the single- 
particle orbitals m\ = S ]X z 2 Ai . Each Fock state can 
be labeled either by A. a partition, or by the occupa- 
tion number configuration n(X) = {n;(A),Z = A$,...,0}, 
where ni(X) is the number of times / appears in A. In the 
spinful case, the partition is replaced by a dressed par- 
tition that mixes momentum and spin. A dressed parti- 
tion (A, a) is given by N e entries Xi and a spin dressing 
o~ € [1, 2, ... C] which obey Aj > Aj+i or A^ = A,+i and 
o~i > cTj+i- The number of groundstates and quasiholes 
states for a given value of N e and iV$ is given by the 
number of dressed partitions (k, r)c admissible. Such a 
partition obeys the following conditions: 



K - K+k > r or Xi - X l+k = r - 1 and <j t < a 



i+k 



(3) 



One can analytically show that for a fixed value of N e 
and JV$, each (k, r)c spinful admissible partition can be 
mapped onto a spinless (fc, C(r— l)+l)i admissible parti- 
tion with the same number of particles and CiV$ orbitals. 
If {(Xi, (Ti)} is a (fc, r)c SU(C) admissible partition, then 
the corresponding spinless partition is {\ = CXi + cr,}. 
In our case, r = 2. Using the relation i'fqh = Cv, we 
obtain that N x .N y = CN$. It is then clear that the 
counting of groundstates and excitations in the energy 
spectrum alone cannot differentiate between spinless and 
spinful (color) FCI states. In particular, a feature of 
these FCI models at fc = 1 has not yet being pointed out 
in previous studies : a (C + l)-fold almost degenerate 
ground state manifold appears at v = 1/(C + 1) for ev- 
ery N x N y values while in the usual spinful FQH (without 
dislocation [45]) one expects to observe them only when 
this number is a multiple of C. This property, which also 
holds true for quasiholes and any fc, would argue for a 
spinless states. 
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FIG. 2. (1,2)3 admissible partitions in the SU(3) case (using 
the red, green and blue colors for each possible value of the 
internal degree of freedom) and their (l,4)i counterparts (in 
gray). Here we have chosen red to be a = 0, green a = 1 and 
blue a = 2. (a) displays a typical admissible configuration for 
the groundstate of the 5(7(3) Halperin state, (b) is an ad- 
missible configuration for the quasihole states that is present 
when computing the PES for the 7V e = 6 particle groundstate. 
(c) is also an admissible configuration for the quasihole states 
but cannot be accessed through the PES from the N e — 6 
particle groundstate. 
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FIG. 3. PES for low energy groundstate manifold on the 
pyrochlore lattice with C = 3 and two body interaction for 
N e = 6 bosons and Na = 3. The number of states below 
the dotted line is 680. The 48 missing states from the (1, 2)3 
counting are the quasihole states that require 3 particles with 
the same value of a. This counting does not match the count- 
ing from the Halperin PES. 



More information about the nature of the ground- 
state can be obtained from the particle entangle- 
ment spectrum [551 [55] (PES). For a d-fold degenerate 
state {\ipi >}, we consider the density matrix p — 
h Si=i iV'i >< V'il- We divide the N e particles into two 
groups A and B with respectively Na and Nb particles. 
Tracing out on the particles that belong to B, we com- 
pute the reduced density matrix pa — Tr^p. The sym- 
metries of the original state, preserved by the operation, 
allows to label the eigenvalues exp(— £) of pa by their 
corresponding quantum numbers. For example in the 
case of a spinful SU(C) eigenstate, both the eigenvalues 
of the Cartan subalgebra and Casimir operators and the 
momentum could be used to label the PES eigenvalues. 
For spinless FQH model states, the number of non zero 
eigenvalues in pa exactly matches the number of quasi- 
hole states for Na < N/2 particles and the same number 
of flux quanta as the original state. The quasihole count - 
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FIG. 4. PES for low energy groundstate manifold on the 
pyrochlore lattice with C = 3 and two body interaction for 
N e = 9 bosons and Na ~ 4. The number of states below 
the dotted line is 14364. This is 1575 states less than the 
(1, 4)i counting. This counting matches the counting obtained 
through the Halperin PES. 



ing is characteristic of each topological spinless state and 
thus the PES acts as a fingerprint of the phase (obtained 
only from the groundstate), also able to differentiate it 
from a CDW, which the energy spectrum is not [34]. In 
FCI, since the groundstates deviate from model wave- 
functions, one expects to observe a low energy structure 
similar to the one of the model state with a gap to higher 
energy excitations. Such a feature has been shown for 
Laughlin and MR- like states in FCI [5] 137], 

When the groundstate possesses an additional symme- 
try preserved by particle partitioning, like the SU (C) 
symmetry for Halperin state in Eq. [T] an additional con- 
straint can reduce the number of eigenvalues of the PES 
from the number of quasiholes of Na particles in the 
original number of fluxes. For instance, Halperin ground- 
states are spin-singlets and hence the number of particles 
per color in them is the same and equal to N e /C. If 
Na > N e /C, the quasihole states with Na particles in a 
given color should not and cannot be found in the FQH 
PES (as depicted in Fig. [2]). In that case, these states 
and their full corresponding SU (C) multiplet have to be 
removed to obtain the PES counting. From this perspec- 
tive, the C = 3 is a perfect test case since C = 2 does 
not provide any additional constraint compared to the 
spinless case due to the constraint on the entanglement 
spectrum Na < N e /2. Indeed, in all the C = 2 cases we 
have checked, the PES counting below the entanglement 
gap matches the spinful (fc, 2) 2 quasihole counting which 
is identical to the spinless (k, 2 + k)\ quasihole counting. 

We now focus on C = 3, k = 1. Figs. [3] and [4] show the 
PES for resp. N e = 6, N A = 3 and N e = 9, N A = 4. 
These system sizes are directly related to their FQH 
counterpart. In both cases, there is a clear gap with 
a counting below it that is lower than the spinless (l,4)i 
counting: this clearly rules out a Laughlin- like v = 1/4 
state. Only the N e = 9 counting matches the Halperin 
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PES counting. The one for N e — 6 is higher than ex- 
pected: it corresponds to the quasihole counting of the 
Halperin state where only the states with more than Na 
particles of the same color are dropped and not the corre- 
sponding full multiplet[2"8]. When N e is not a multiple of 
3, the PES exhibits an entanglement gap above a (l,4)i 
counting. 

In conclusion, we presented strong evidence for the ex- 
istence of a series of abelian and non-abelian states at 
filling v — k/(C + 1) in fractional Chern insulators with 
single bands of Chern number C > 1. The energy spectra 
for both ground-state and quasiholes suggest that several 
models support these clustered states. In order to iden- 
tify the states, we presented evidence from the entan- 
glement spectrum of a hidden SU (C) symmetry of the 
FCI states that seems to rule out a spinless description. 
However, we observe several rate anomalies (at particu- 
lar commensurate particle numbers) in the FCI entan- 
glement spectrum compared to that obtained from the 
FQH, which could be due to the presence of "effective" 
twisted boundary conditions in the FCI problem. 
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